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Abstract 



Using an isomorphism between Hilbert spaces I? and I 2 we consider Hamiltonians 
which have tridiagonal matrix representations (Jacobi matrices) in a discrete basis and 
an eigenvalue problem is reduced to solving a three term difference equation. Tech- 
nique of intertwining operators is applied to creating new families of exactly solvable 
Jacobi matrices. It is shown that any thus obtained Jacobi matrix gives rise to a new 
exactly solvable non-local potential of the Schrodinger equation. We also show that 



the algebraic structure underlying our approach corresponds to supersymmetry. Su- 



percharge operators acting in the space £ 2 x £ 2 are introduced which together with a 
matrix form of the superhamiltonian close the simplest superalgebra. 



1 Introduction 



Supersymmetric quantum mechanics (SUSY QM) first introduced by Witten pQ as a simplest 
super symmetric model of the quantum field theory represents now a fast progressing field of 
modern theoretical physics (see recent books 0). As it has been understood after the work 
by Andrianov et al [3], the SUSY QM may be considered as another formulation of Darboux 
transformation method well-known in mathematics from the original paper by Darboux 0] , 
book by Ince {5^ , and book by Matveev and Salle [U] where this method is widely used in the 
context of the soliton theory. An essential ingredient of the method is a particular choice of 
a transformation operator (or intertwining operator [7j) in the form of a differential operator 
(see e.g. jS]) which intertwines two Hamiltonians and relates their eigenfunctions. This 
approach being applied in quantum theory allows one to generate a huge family of exactly 
solvable local potentials starting with a given exactly solvable local potential. 
Published in Phys. Lett. A 302 (2002) 234-241 
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The technique of transformation operators, first proposed by Delsart [5| (he calls them 
transmutation operators, see also e.g. [?]) in connection with generalized translation oper- 
ators, is essentially based on intertwining relations, which have an universal character [TUj . 
It is applicable not only to differential equations (like the Schrodinger equation) but also to 
difference equations |TT|. IT2*1 |H] which, in the simplest case, represent eigenvalue problems for 
Jacobi matrices. Intertwining operators are widely used in the context of soliton solutions 
of nonlinear lattice equations (see [S] and references therein). 

It is known that spectral problems for finite-difference (or discrete) Schrodinger equations 
have much in common with the theory of orthogonal polynomials The Darboux trans- 
formation method gives new insight to this theory and is very fruitful in establishing new 
properties and deeper understanding some known relations [T3|. As it was shown in the 
continuous Darboux transformation is intimately related with quantum inverse scattering 
problem. A discrete version of this problem has been considered in (for a review 

see e.g, jTH]) but its relationship with discrete Darboux transformations is unknown to the 
authors. We note that considerable attention has been drown to the non-Hermitian form 
of the problem @ [TTJ, [T2J [Hj, while in quantum mechanics one usually uses self-adjoint 
Hamiltonians. The study of the Hermitian problem has been just started [IZ|. We thus 
notice that the discrete version of the Darboux transform though ascending to Christoffel 
[TT] and Geronimus works is developed much less than the continuous one. 

This paper is aimed to partially fill in this gap. We apply the intertwining relation to a 
discrete eigenvalue problem for a Jacobi matrix. This allows us to generate a family of new 
exactly solvable Jacobi matrices. Eigenvalue problems for these matrices naturally appear 
in quantum mechanics when a discrete basis is used for solving the Schrodinger equation. 
In particular, there exist exactly solvable Hamiltonians represented by Jacobi matrices [T5] . 
We show that new Jacobi matrices define new exactly solvable potentials of a non-local 
nature. Moreover, when two lattice eigenvalue problems are written as a single matrix 
problem, nilpotent difference supercharges may be introduced to map the solutions of these 
problems to each other. We also show that these supercharges together with a difference 
superhamiltonian close the simplest superalgebra. We associate these supercharges with 
discrete supersymmetries of the Schrodinger equation. 

2 Intertwining in ^ 2 -space 

Consider self-adjoint Hamiltonians ho and hi defined in a Hilbert space H = L 2 (IR). We do 
not exclude the presence of a continuum spectrum and if necessary consider them as defined 
in a wider space H_ of linear functionals over H + , H + C H C f/_ , where H_ and H + are the 
Hilbert-Schmidt equipment of the space H (so called Gel'fand triplet, see e.g. [THj) without 
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especially mentioning it. Let \iPe) and \iPe) be eigenkets of ho and hi respectively with an 
eigenvalue E, h \ip E ) = E\ip E ), hi\ip E ) = E\ip E ), and \n) be a basis (orthogonal or not) in 
H. The basis \n) is supposed to be such that the action of both ho and hi takes the form of 
three term relations: 

h \n) = a n \n - 1) + a n+ i\n + 1) + q n \n) (1) 

hi \n) = a n \n - 1) + a n+ i\n + 1) + q n \n) . (2) 

We assume that n takes positive integers including zero, n = 0, 1, 2, . . . and a = a = 0. 

Denote by ?/>„ = ip n (E) and ■?/>„ = ip n (E) the Fourier coefficients of \ip E ) and m terms 
of the basis \n) respectively: 

\^e) = ^ip n \n) , fe) = 5^^n|n). (3) 

Equations (0) and © define the action of ho and /ii on any = Yl c n\ n ) ^ H where the 
sum contains a finite number of items. (They form a finite set in H .) An eigenket of a 
Hamiltonian belongs usually to a wider domain which can be specified in terms of its Fourier 
coefficients over a basis. We will not discuss these mathematical subtleties here and shall 
simply suppose that all Fourier series converge and when an operator acts on a series it may 
be moved through the sign of sum. Thus, having (0) and (j2J in mind we can act by h and 
hi on \ip E ) and \ip E ) given by Q respectively. As a result we obtain the Fourier series 

W\^e) = ^2(h ip) n \n) , hltpE) = y^(fei-0) n |n) (4) 

n n 

with the Fourier coefficients 

(h 4>) n = a n Vn-i + a n+ itp n+ i + q n tp n , (5) 
{hii)) n = a n 4) n _i + a n+1 ^ n+1 + q n i) n . (6) 

Note that these quantities depend on the energy E but we do not indicate this to avoid 
cumbersome notations. The eigenvalue problems for h and hi are reduced now to the 
second order three-term finite difference equations 

a n ip n _i + a n+ iip n+ i + q n i) n = Eip n , (7) 
a n ?/v_i + a n+ iV>„+i + Qni'n = Ei> n . (8) 

Note that the space of sequences {ip n } with an appropriately defined inner product is 
known as the space £ 2 which also is a Hilbert space (if the space is not complete, it always 
may be enlarged to become a Hilbert space). If necessary this space may also be equipped 
with t±, t\ C £ 2 C l\ The meaning of £± is the same as H±. Eqs. © and © define 
matrix representations of the Hamiltonians ho and hi. 
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A second order finite difference equation like a second order differential equation has two 
linearly independent solutions. This means that for a fixed E one has a two-dimensional 
linear space of solutions to Eq.([ZI) (and respectively (jHJ))- When the spectrum of a Hamil- 
tonian is non-degenerate and the energy is fixed, in this space may exist a unique (up to 
normalization) element that can be associated with an eigenvector of ho- For a bound state 
eigenvalue E it may be selected from the condition 



(h 4> E \hoi)E) = ^2(h ^)*k (ho^) n (k\n) < oo . 



n, k 



We use "*" to distinguish complex conjugate quantities. For a continuous spectrum this 
value may diverge like the Dirac 5 function, 

(MfcMis') = X>o^)fc(£) • (M0n(£')(%> ~ 8{E — E') . 

n,k 

In these equations (hoip) n = (h Q ip) n (E) are defined in (jSJ). Note also that the continuous 
spectrum may be two-fold degenerate. In this case any solution of Eq. ((7j) has a physical 
meaning. The right-hand sides of Eqs. (JSJ) and © define the action of the operators ho and 
h\ in the space I 2 , e.g. h {ip n } = {(h ip) n }. 

Denote ip n another solution of Eq. (J7|) corresponding to the same energy E. Then 
excluding q n from eigenvalue equations for ijj n and ip n , one arrives at the equation 

OnO/vAn-l - fn-lVVi) = Qn+lft+lV'n ~ Vv/Wl) (9) 

which gives us a discrete analogue of the Wronskian as a combination of ip n and ip n inde- 
pendent of the discrete variable n 

W^n, 1p n ) = On^Vn-l ~ ^n-l^n) = W = COnst(n) . (10) 

When ivq and ip n are fixed, this equation is a recursion relation for i[) n which has the following 
solution ^ 

For solutions ip n and £ n of (J7J) with different eigenvalues E and A respectively similar calcu- 
lation gives a recursion relation for W(£ n ,i/j n ) = W n , 

W n+1 = W n + (X-E)^ n , (12) 

which may be easily iterated to become 

n 

W n+ i = W l + (X-E)J2 ikipu ■ (13) 

fc=i 
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Let us suppose that the eigenkets of ho, \?Pe), are known, i.e. one knows solutions of Eq. 
(J7J). Then the eigenvalue problem for hi, i.e. the search for solutions of (jHJ), may be replaced 
by the search for an intertwiner L defined as 

Lh = hiL. (14) 

Once L is found, the eigenkets \iJje) are obtained by applying L to \4>e), \^e) = L\^e)- 

Let us consider a particular anzats for L assuming that it is defined in terms of the basis 
|n) as follows: 

L\n) = A n \n- 1} +B n \n) . (15) 

Then the Fourier coefficients ip n = Lip n = (Li[)) n of \ip) = L\ip) = X^n^nl^) are expressed in 
terms of i]) n : 

%jj n = (Lip) n = A n+1 ip n+1 + B n ip n . (16) 

To determine A n , B n defining L and a n , q n defining hi, we shall use the intertwining 
relation ([T3)l. The successive application of © and (fTB^) gives 

L{h ip) n = A n+1 (a n+2 ip n+2 + q n+ i^n+i + a n+1 ^ n ) + B n (a n+1 ip n+1 + q n ip n + a„/0n-i) 

whereas the formulae (fTB^l and © yields 

from which one finds the following system of equations: 

A n o. n +i — A n+ ia n , (17) 

B n a n = B n _ia n , (18) 

A n +iq n +i + B n a n+ \ = B n+ ia n+ i + A n+ iq n , (19) 

A n +ia n+ i + B n q n = a n A n + B n q n . (20) 

Having found a n from (fT8|) and replaced it in (JT7J), one can "integrate" the latter equation 
to get 

B n = A^ (21) 

where A denotes an "integration constant". After excluding q n , a n and B n from (|19|l and 
(|2()j) one obtains the equation for A n 



Aa n+ i A n _ Aa n+ 2 A n+i 

Qn 72 a n+l "j 0,71 ~ Q n+1 72 a n+2 ~7 O-n+1 ■ \**) 

A n+1 Aa n A n+2 Aa n+ \ 



Clearly, both left- and right-hand sides of this equation are independent on n. Then denoting 
this quantity by A and introducing a new variable £ n 

£n-l A n 



£n Aa n 
5 



(23) 



one arrives at the equation for £ n : 

£n-lOn + + (q n - A)£ n = . (24) 

The latter equation is nothing but the initial eigenvalue problem (j7J) which is supposed to 
be solved. We thus shall suppose that the sequences {£ n (A)} and {£n(A)} being linearly- 
independent solutions of (|24jl for a fixed value of A are known. When A = E is a point 
of the spectrum of h , there exists their linear combination which gives the ket \iPe) in n- 
representation for non-degenerate E. For two-fold degenerate E both {£ n (E)} and {£ n (E) } 
define two linearly independent kets. 
Having found A n and B n : 

A n = {-Aa n i n ^li n \ 1 / 2 , B n = -\-Aa n+1 i n+1 /i n } 1 ' 2 , (25) 

one gets solutions ip n (E) of the transformed equation (JSJ) 

-Aa n+ i 



ME) = (Lip) n (E) 



1/2 



{^ n+ i{E) - in+iME)) , (26) 



£n£n+l 

{h^) n {E) = Ei, n {E) . (27) 
From (|18j) one obtains the elements a n 

0-n = T - [Ondn+l^n+l^n-l] 1 ^ 2 (28) 

and from (|19|) and (|2U|) one has two equivalent expressions for q n : 

~ | Sn— 1 £n _ £n+l £ra+2 fnrA 

Qn — In T Ctn - T «n+17 — <?n+l — a n+l — T h «n+27 • (/^J 

The formula ()26j) defines the sequence {^ n (i2)} for any E ^ X. When E = A, it gives 
zero, {L^) n = 0, but we can find a solution of Eq. (JHJ) at £7 = A by acting with L to another 
solution of Eq. (|24j). £„, which may be found using (fTTJ) with the replacement ip — ► £. Thus 
we have the sequence given by 

r? n = (L£)„ = [a n+ i£ n £ n+ i]~ 1/2 , (hir]) n = Xr] n (30) 

where inessential constant w is omitted. Another solution rj n of Eq. (JHJ) at = A may be 
found by appropriately changed formula (fTTj) 

n 

rjn = Vn [Vo/Vo + w ^2^ k ] , (h 1 r}) n = \rj n . (31) 
k=i 

So, we know all solutions of the transformed discrete equation (jHJ) with a n and q n given by 
(I28|) and (J29j) provided all solutions of (J7J) are known. When E ^ X, they are given by 
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(J2fijl and for E = A the formulae (J3Uj) and (j3~Tj) should be used. Note that the fact that the 
sequence {ip n } = {{Lip) n } satisfies the difference equation (jHJ) means that in the space £ 2 
the intertwining relation Lho = h\L holds where the action of ho on a sequence is defined 
by (0) and the action of L is defined by (fH)|) with A n and B n from 

The most popular bases in quantum mechanics are orthonormal bases, {n\n') = 8 n ,n'- 
For such a basis the action of the operator L + , Hermitian conjugate to L with respect to 
the inner product, on any finite vector \ip ) = Yl n ^nl 71 ) where the sum contains only a finite 
number of items can be easily found. For this purpose one first notices that for an L given by 
(fTKj) one has (n'\Ln) = A n 5 n > tn -i + B n 5 n > tn = A n i + iS n ' + i jn + B n 5 n r, n . Now from the condition 
(n'\Ln) = (L + n'\n) one deduces that L + \n) = A^ +l \n + 1) + B^\n) from which it follows 
that 

L + $) =J2( L+ ^)n\n} } (L + ^) n = A*J n - 1 + B*J n (32) 

n 

with A n and B n given in ()25|). Since the Hamiltonians h and hi are self-adjoint, the 
conjugate operator L + satisfies the conjugate intertwining relation h$L + = L + h\ which 
means that L + realizes the transformation in the opposite direction, from solutions of (JBJ) to 
solutions of (|7j) but it is not inverse to L. We shall find the superpositions L + L and LL + in 
the next section. 

For a self-adjoint Hamiltonian it is natural to suppose that a n and q n are real for all 
n = 0, 1, . . .. This implies that for a real A all £ n and £ n may be chosen real which from 
now on will be assumed. This means that the variables A n and B n defined by ()25j) may 
become complex only owing to the constant —A. Since the operator L + is well-defined for 
any finite element from £ 2 , it is not difficult to convince ourselves that (L + rj) n = and 
(L + ff) n = —woV—A* £ n . Since the usual condition for a Wronskian of two solutions 
belonging to the same energy is Wq — 1, the natural choice for A is A = — 1 . In this case one 
recovers a discrete analogue of the continuous Darboux transform (L + rj) n = — £ n . 

So, we see that the discrete Darboux operators L and L + realize a one-to-one mapping 
of the spaces of solutions of discrete Schrodinger equations (0 and (jHJ) for any E ^ A. 
Moreover, this mapping may be continued to E = A by putting £ n < — > rj n and £ n < — > r\ n and 
considering this as a linear mapping. We conclude hence that the knowledge of all solutions of 
the initial equation (J7|) provides the knowledge of all solutions of the transformed equation (jHJ) 
and, in particular, one can get all "physical" solutions in the usual quantum mechanical sense. 
For this purpose one notices first that the norm of the vector ip = {ip n } = Lip = {Lip n } is 
proportional to the norm of the vector ip = {ip n }, (ip\ip) = Ylm 1 Pn l l J n = ^ n C(^ + ^)n = 
(E — A) J2 n ^n^Pn = (E — A) where the use of the factorization property has been made. 
This means that all "physical" solutions of (JTj) with E ^ A are mapped by L onto "physical" 
solutions of (jHJ). It is evident that the inverse statement is also true but instead of L one has 
to use L + . To find all "physical" solutions of (jHJ), it remains hence to analyze the sequences 
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with eigenvalue E = X forming a 2- dimensional space spanned by {rj n } and {rj n } which is 
not a difficult problem. 

As a final remark of this section we note that the formulae (fT6|) and f!32[) define operators 
L and L + only on finite elements from H (finite linear combinations of basis vectors \n)). But 
using the technique of polar decompositions of closed operators in terms of quasi-projectors 
developed in [20] for a continuous basis, one can extend this action such that the operators L 
and L + become mutually conjugated with respect to the inner product in H with well-defined 
domains of definitions. 



3 Discrete supersymmetries of the Schrodinger 
equation 

The usual supersymmetry of the continuous Schrodinger equation is based on factorization 
properties of Darboux transformation operators (see e.g. jHj). Let us find their counterparts 
for the discrete equations (0) and ©• For this purpose, let us consider the action of the 
superposition of L and L + on a ket- vector Using (fTp^) and (J3~2")l one finds 

L + L\^)=Y J {L + L^) n \n), 

n 

{L + Li,) n = A* n B n ^ n ^ + B*A n+1 iP n+1 + (\A n \ 2 + \B n \ 2 )^ n (33) 

where tp n are Fourier coefficients of a vector e H over the basis \n). If for the sake 
of definiteness we assume that A < q n for all n, then from Eq. h follows that both 

dnin-i/in < and a„+i£„+i/£ n < which permits us evaluate absolute values in (J33J). 
Finally, using once again the equation for u n one g e ts 

(L + L^) n = \A\ ( [ho - XWn ■ (34) 

Similar calculations lead to another factorization 

(LL + ip) n = \A\([h 1 -X}^) n . (35) 

It is necessary to note that these factorizations take place for any sequences {ip n } and {ip n } 
and the superposition of operators L and L + acting in the space of sequences £ 2 is defined in 
the regular way L + Lip n = L + (Lip) n = L + ip n = (L + ip) n . But when {ipn} is an eigenelement of 
ho, hotp n = (hoip) n = Eip n , one gets from (pMjl (L + Lip) n — (E — X)ijj n where we put A = —1. 
Similarly, when hiip n = (h\i})) n = Eip n one obtains using (|33|) (LL + ip ) n = (E — X)^ n 

We can now imitate supercharge operators by introducing nilpotent difference matrix 
operators 



S 



acting in a vector space of 2-component column-vectors 



where ip n and ip n are elements of arbitrary sequences from i 2 . The superhamiltonian 

together with supercharges (}3T)j) close a simplest superalgebra 

[Q,H]= [Q+H] =0, Q 2 = (Q+) 2 = 0, {Q,Q+}=H-AI (39) 

where I is the 2x2 unity matrix. Here commutation relations follow from the intertwining 
relation for L, L + , ho and hi and the anticommutation relation is an implication of the 
factorization properties for these operators. The discrete eigenvalue problem for the super- 
hamiltonian HI is nothing but the collection of either equations (J7J) and © in a single matrix 
equation. Since the operators Q and Q + from (|3*9*j) act in the space £ 2 , we associate them 
with discrete supersymmetries of the Schrodinger equation. 



4 New non-local exactly solvable potentials 

Now we shall apply the technique of discrete SUSY (or equivalently discrete Darboux trans- 
formations) to generating non-local exactly solvable potentials. 

Let us consider the free particle Hamiltonian h Q = p 2 = —d 2 /dx 2 . Since the momentum 
operator may be expressed in terms of the harmonic oscillator creation a + = id/dx + ix/2 
and annihilation a = id/dx — ix/2 operators p x = —(a + a + )/2, the Hamiltonian ho is a 
quadratic form of a and a + , ho = (a + a + ) 2 /4. Therefore, the action of ho on the oscillator 
basis \n) which in coordinate representation looks like 

V>»(aO = (x\n) = {-i) n {n\2 n V2^)- l ' 2 e- x2/i H n {x/V2) (40) 

takes the form of a three term relation 

ho\n) = yn{n-l)\n - 2) + \ ^{n + l)(n + 2)|n + 2) + (f + \)\n) , (41) 

where the use of the well-known properties of the creation and annihilation operators a\n) = 
\fn\n — 1), a + \n) = \Jn + l|n + 1) has been made. Let \ipE) be a continuous spectrum 
eigenket of ho, ^oIV'b) = E\iPe)- Then using the self-adjointness of ho we get for the inner 
product ip n = ip n (E) = (ipE\n) the following discrete eigenvalue problem 

a n 1p n - 2 + a n+2 1p n +2 + Qn^n = Eip n , (42) 

a n = yn(n-l), 9n =(§ + 1). (43) 
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A "physical" solution to this problem, ip n = ip n (E), can be easily obtained since it coincides 
with the Fourier image of the function (|40p 

A = 2(n\2 n V^)- 1/2 e~ E H n (V2E) . (44) 



Here H n (z) are Hermite polynomials. It is not difficult to convince ourselves that Eq 
represents another form of the well-known recursion relations for the Hermite polynomials. 

The Hamiltonian h is the kinetic energy operator for h = hi = h + V . Let us consider 
V such that 

V\n) = d n \n - 2) + d n+2 \n + 2) + r n \n) (45) 

In coordinate representation this operator is not reduced to the multiplication on a function. 
We conclude therefore that this is a non-local operator. 

Let \4>e) be an eigenvector of h, h\ipE) — E\i/je)- Then the eigenvalue problem for the 
inner products ip n = i/j n (E) = (i/ieIA) is just reduced to the discrete eigenvalue equation 
(jHJ) which in our case has the form: 

a n i> n - 2 + a n+2 ^ n+ 2 + q n ^n = Etjj n (46) 

where 



a n = d n + \\Jn{n - 1) , q n = r n + § + \ . (47) 

Consider now a subclass of potentials (J45|) for which Eq. (|46|) coincides with the Darboux 
transform of the Eq. (|4*2*|) . In this case the functions a n and q n are defined by Eqs. (|2*H|l and 



«n — [(ln a n+2^n-2^n+2\ ^V^n, (48) 
Qn = Qn+2 — «n+2^n+2/ £n + 0"n+i^,n+4:/ £,n+2 (49) 



where a n and q n are given in (|45|) and ^ n is a particular solution of (j42j) at E = A < 0, i.e. 

i n = {n\2 n y l l 2 H n {V2X) . (50) 

Note that these functions are real for n even and purely imaginary for n odd and therefore 
a n and q n are real in either case. Solutions of Eq. (|46J) are found with the aid of Eq. 



Q"n+2 



n+2 



1/2 

(CnA+2 ~ in+2^n) (51) 



which gives us continuous spectrum eigenfunctions for the Hamiltonian h = ho + V with the 
interaction V given in the form of an infinite tridiagonal matrix 

(k\V\n) = d n 6 k>n - 2 + d n+2 5 k>n+2 + r n 5 k>n (52) 
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with d n and r n denned with the aid of (j4Tj) . (jlHj) and (jUJ). It is not difficult to see that r\ n 
found by (j3*Uj) behaves as n~ l l 2 and given by (|31|) behaves as n 3 / 2 at large n. We conclude 
hence that no discrete levels in the potential V is created by our method and this is a purely 
scattering potential. 
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